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ABSTRACT 

Fermi has resolved several star-forming galaxies, but the vast majority of 
the star-forming universe is unresolved and thus contributes to the extragalactic 
gamma ray background (EGB). Here, we calculate the contribution from star- 
forming galaxies to the EGB in the Fermi range from 100 MeV to 100 GeV, due 
to inverse-Compton (IC) scattering of the interstellar photon field by cosmic-ray 
electrons. We first construct a one-zone model for a single star-forming galaxy, 
assuming supernovae power the acceleration of cosmic rays. The same IC inter- 
actions leading to gamma rays also substantially contribute to the energy loss of 
the high-energy cosmic-ray electrons. Consequently, a galaxy's IC emission is de- 
termined by the relative importance of IC losses in the cosmic-ray electron energy 
budget ("partial calorimetry" ) . We use our template for galactic IC luminosity 
to find the cosmo logical contribution of star- forming galaxies to the EGB. For all 
of our models, we find the IC EGB contribution is almost an order of magnitude 
less than the peak of the emission due to cosmic-ray ion interactions (mostly 
pionic PcrPism — > 7r° — > 77); even at the highest Fermi energies, IC is subdomi- 
nant. Moreover, the flatter IC spectrum increases the high-energy signal of the 
pionic+IC sum, bringing it into better agreement with the EGB spectral index 
observed by Fermi . Partial calorimetry ensures that the overall IC signal is well 
constrained, with only modest uncertainties in the amplitude and spectral shape 
for plausible model choices. Partial calorimetry of cosmic-ray electrons should 
hold true in both normal and starburst galaxies, and thus we include starbursts 
in our calculation. We conclude with a brief discussion on how the pionic spectral 
feature and other methods can be used to measure the star-forming component 
of the EGB. 
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Introduction 



The window to the high-energy (> 30 MeV) ga mma-ray cosmos has be en open now for 
four decades, with measurements by OSO-3 satellite (IKraushaar et al.lll972l ) followed by the 



second Small Astronomy Satellite (S AS-2)(Fichtel et a 



Ray Experiment Telescope (EGRET) (jSreekumar et al 



19751) and the Energetic Gamma 



19981 ). These revealed the existence 



of a diffuse, extra-Galactic gamma ray background (EGB). More recently, the advent of 
the Fermi Gamma-Ray Space Telescope has substantially sharpened our observational view 
of the EGB. With better energy and angular resolution and much higher sensitivity than 
EGRET, Fermi has resolved many more gamma-ray point sources and better determined 
the diffuse background and its energy dependence. The EGB data are consistent with a 
power law of spectral index 2.41 ± 0.05 ( lAbdo et al.ll2010if ) for energies > 100 MeV. 



The origin of the EGB remains an open question. A contribution from active galaxies 



(e.g., 


Strong et al. 


1976a 


Abdo et al. 


2009: 


Stecker & Venters! 


201 




\bazaiian et al. 


2011) 


and star- forming galaxies (e.g., 


Strong et a . 


1976b 


Bignami et al. 


1979: 


Pavlidou & Fields 


2002; 


Fields et al. 


2010; 


Stecker & Venters 


2011 




Vlakiva et al. 


2 


01 


.) is "guaranteed" in 



the sense that these are known, resolved extragalactic source classes that must have un- 



resolved counterparts that will contribute to the EGB. Oth er possible 



truly diffuse emission such as dark matter annihilation (lAndo et al. 



from cosmic rays accelerated in structure formation shock s (e.g. 



Miniatil 12002). unresolved ordinary and millisecond pulsars ( IFaucher-Giguere &: Loeb 



Geringer-Sameth &: Koushiappas 



1GB sources include 



20071 ) , intera ctions 
2~00~ol : 



Loeb &: Waxman 



20 



interactions with Oort cloud bodies ( IMoskalenko fc Porter! 120091 ) 



20101 : 



2 ), and even Solar-System emission from cosmic-ray 



One of Fermi's major achievements has been to establish external star-forming galaxies 
as a new class of gamma-ray sources. These detections give a global view of the gamma- 
ray output as a result of star-formation, complementary to the resolved but local images of 
the Milky Way. Fe rmi has not only d etected but also spatially resolved the La rge Magel- 



lanic Cloud fLMC; lAbdo et a 



2010eh. As antic ipated flPavlidou & Fields! 120011 ). the SMC 



( lAbdo et al.ll2010b[ ) and M31( ]Abdo et al.ll2010dl ) have also been detecte d, while other nor 



mal star-forming galax ies in the Local Group (including M33) have not ( I Abdo et al.ll2010dl : 
Lenain &; Walter! 1201 if ). Beyond the Local Group, Fermi has detected starbursts galaxies 
characterised by very high star- formation rates, as anticipated bylTorres et al.l ((2004J). Fermi 
has d etected the starbur sts M82 and NGC 253 kbdo et al.lboiOah and NGC 1068 and NGC 
4945 (TNolan et al.ll2012l ). 



The Fermi star-forming galaxies offer a qualitately new probe of cosmic rays; they also 
inform and calibrate efforts such as ours to understand the EGB contribution from the vast 
bulk of the star-forming universe that remains unresolved. The LMC is the best resolved 
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individual system, and there the energy spectrum is consistent with pionic, while t he spa- 
tial distribution can be used to study cosmic-ray propagation (IMurphy et al.ll2012l ). More 
broadly, the ensemble of all Fermi star-forming galaxies encodes information about global 



cosmic-ray energetics an d interaction mechanisms (IPersic fc Rephaefill2010l ; lLacki et al.ll2011 



Persic fc Rephaefill2012l ). In particular, Fermi reveals a strong correlation between gamma- 
ray luminosity L 1 and supernova rate (or equivalently star- formation rate ip). This is ex- 
pected if supernovae provide the engines of cosmic-ray acceleration. Remarkably, all star- 
forming galaxies d etected to date can be well-fit with a single power law L 7 oc if) 



1.4±0.3 



Abdo et al.ll2010dh 



e.g. 



The main mechanism of gamma ray production in star-forming galaxies is anticipated 
to be the same that dominates Milky Way diffuse gamma rays: pionic emission p CT pi sni 
PP7T — > 77, arising from interactions between cosmic-ray hadrons (ions) and inte rstellar 
gas Jstecker & Ventersl lioiH I Abdo etaDbood : [Fields et al.l lioiol : Istrong et al.lboioh . This 
mechanism is likely responsible for the non-linear relation between the luminosity of Fermi 
galaxies a nd their star-formation rate. Namely , the observed correlation is co nsistent with 
a picture JPavlidou k Fieldsl l200ll : iFields et aliboid : IPersic k Rephaelill201lh in which the 
cosmic-ray proton flux is controlled by the supernova rate, the total number of targets is set 
by the ga laxy's gas mass, and the gas and star-formation are linked by the Schmidt-Kennicutt 
relation (jKennicuttlll998l ). 



Here, we examine the cosmological contribution of star-forming galaxies by the inverse- 
Compton (IC) scattering e cr 7i S — > cy of cosmic-ray electrons on the interstellar radiation 
field (ISRF). In order to do so, we construct a one-zone model of a star- forming galaxy, 
with the IC emissio n normalised to the Milky Way IC emission as computed in GALPROP 
(IStrong et al.ll2010l ) and use this as our galaxy template. 



The IC gamma r ays are produced by upscattering of interstellar radiation by high -energy 
cosmic-ray electrons ( Felten Sz Morriso"nlll963 ; Felten 1965 ; Brecher fc Morrison 1967 ). Inverse- 
Compton scattering also represents an important energy loss mechanism for these electrons; 
the other important losses are bremsstrahlung and synchrotron. The relative importance 
of these losses depends on the cosmic-ray energy and on interstellar radiation and matter 
densities. Where inverse Compton losses dominate, the energy injected into cosmic-ray elec- 
trons is ultimately re-emitted as IC gamma-ray photons. This equality of energy loss and 
energy outp ut is known a s calorimetry, and was first described in the context of high-energy 
electrons bv lFeltenl (119651) and explored in detail for diffuse high-energy gamma rays by cos- 
mic rays by iPohll (11994I ). If the other losses are negligible then we have perfect calorimetry. 
However, if other losses compete, the gamma-ray energy output is reduced by the IC fraction 
of the energy losses by the cosmic rays. As we will see, the latter holds true for the case of 
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Milky Way type galaxies. We can explain our results in terms of this partial calorimetry. 

In contrast, cosmic-ray hadrons (mainly protons, as well as other ions) in the Milky 
Way suffer losses dominated by escape rather than collisions. The pionic emission from 
normal star-forming galaxies is thus not calorimetric; however, in starburst galaxies , proton 



inelastic interactions can dominate losses and lead to calorimetry ( jLacki et al .1120 101 ). In the 



Milky Way, pionic emission dominates the total Galactic gamma-ray output (luminosity) , 



excee ding IC emission by factors of up to ~ 5 in GALPROP calculations (IStrong et al. 



20101 ). Our analysis shows that the inverse-Compton component from star- forming galaxies 
over cosmological volumes is nearly an order of magnitude lower than the peak of the pionic 
component. 

The paper is organised as follows. Section [2] gives an order-of-magnitude estimate of the 
inverse Compton EGB contribution by star-forming galaxies, and serves as an overview to 
our paper. We first build a template for the IC emission from a star-forming galaxy, starting 
with the various components of the background interstellar photon field These interstel- 
lar photons serve as scattering targets for the cosmic ray electrons, whose propagation is 
explained in §H For the highest energy background photons and electrons, the scattering 
occurs in the Klein-Nishina regime, which affects the emergent spectrum and is detailed 
in §5.11 Our one-zone model for the inverse-Compton luminosity from a single galaxy is 
presented in §5.21 In £j6l the total intensity over cosmological volumes is calculated and 
compared with the pionic component. Section [7] discusses the implications of our results. 



2. Order-of-Magnitude Expectations 

An order-of-magnitude estimate of our final result will help frame key physical issues 
and astrophysical inputs. Our goal is to find the gamma-ray specific intensity, Ie, due to 
inverse Compton scattering in star-forming galaxies, at energy E. For photons up to at 
least < 30 GeV, the universe is optically thin; thus, the intensity is simply given by the 
line-of-sight integral 



where £ 7 is the luminosity density or cosmic volume emissivity of the galaxies, and g?h — 
c/Hq = 3000/i -1 Mpc is the Hubble length. The total luminosity density of a distribution of 
IC-emitting, star-forming galaxies can be expressed as a product of the luminosity function, 
drigax/dL^ and the luminosity of each individual galaxy, L 1 : 



(2) 
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where Le is the luminosity of an average galaxy at energy E. 

Within a single galaxy, the gamma-ray luminosity is an integral 



L~f(E e 



<?ic dV e 



(3) 



over the volume in which cosmic-ray electrons propagate. Here the inverse- Compton volume 
emissivity q ic depends on the product of the targets and projectile densities, and the inter- 
action cross section. The density of targets is simply the number density of the interstellar 
photons, rii sr f. The projectiles are cosmic-ray electrons, with flux density <fi c . The cross sec- 
tion is that for inverse Compton scattering, da- ic / dE em . Therefore, the emissivity for a single 
galaxy is given as, 



'dK, 



-n isr f 



(4) 



with the angle brackets indicating averaging over the cosmic-ray and background photon 
energies. 

The re is substantial evidence that cosmic-ray acceleration is powered by s upernova ex- 
plosions (IHayakawa et al.lll958l . Il964i lUchiyama et al.ll2007t lAhlers et al.ll2009l ). Indeed, the 



cosmic-ray / supernova lin 



via radio (Webber et al. 



t historically has been b e tter established for cos mic ray electrons 



1980l: Istrong et al. 



chrotron flUchivama et al 



Fermi -LAT (lAhlers et al 



2007: 



Helder et al 



2011 



Bringmann et al.l 120121 ) and X-ray syn 



20091 ) measurements. Only very recently have 



2009 



I) measurements of pionic gamma rays provided the most di- 



rect evidence for supernova acceleration of protons and other ions. Consequently, a galaxy's 
cosmic-ray injection rate is proportional to its supernova rate, q cr oc i?sN and therefore, the 
galactic star-formation rate, if). 

Cosmic-ray electron propagation is dominated by energy losses in the form of inverse 
Compton, bremsstrahlung and synchrotron processes. Each contributes to a total energy 
loss rate b tot = \dE e /dt\. The IC loss rate is proportional to the background photon energy 
density: b ic oc U isT { oc n isr f. In steady-state, cosmic-ray losses balance their production, 
and thus the flux is set by the production rate q e times the stopping time oc l/b tot , giving 
(fi e oc q e /btot- If IC losses dominate, then (fi e oc q e /b\ c oc g e /^isrf, and the cosmic-ray flux 
is inversely proportional to the photon density: more interstellar photons mean a shorter 
stopping time. More generally, the cosmic-ray flux (fi e oc (b lc /b tot ) q e /n isri is lower by the 
fraction b ic /b tot of energy losses in IC. 

The galactic IC emissivity depends on the product of flux and targets: q\ c oc e n isr f oc 
(&ic/frtot)<? e - If IC losses dominate, then 6 tot ~ b- lC) and we arrive at the simple and important 
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result q- lc oc q e : the IC photon production rate is proportional to the electron injection rate. 
Physically, this arises because when IC losses dominate, the cosmic ray losses are directly 
probed by the gamma-ray signal we are calculating. Thus, the conservation of energy relates 
the portion of energy injected into cosmic-ray electrons to that later emitted as IC photons. 
Hence, the IC emission serves as an electron "calorimeter." This limit is the case of perfect 
calorimetry, but in reality synchrotron and bremsstrahlung losses are also present and can 
be non- negligible. This is the case of partial or fractional calorimetry where q\ r J q P . oc b\ r Jb tn t: 
IC photons trace the portion of cosmic-ray energy lost via this mechanism (jPohll Il994l ). 
As we will see in greater detail in §HJ IC losses are always important yet do not vastly 
overwhelm the other losses. Thus, calorimetry is only partially realised in detail. But even 
the approximate validity of calorimetry is sufficient that the IC gamma ray luminosity is a 
fairly robust calculation. 

In the limit of perfect calorimetry, the IC volume emissivity is qi C oc q e , and integration 
over all of the supernovae acting as cosmic-ray accelerators gives the IC luminosity Le oc 
.Rsn oc ip- Physically, a fixed fraction of each supernova's energy (and thus a fraction of each 
parcel of mass in new stars) goes into cosmic-ray electrons and ultimately into IC photons. 
This scalin g can then be calib rated by GALPROP estimates for the total Milky- Way IC 
luminosity fjStrong et al.l 120101 ). which determines the IC output per unit star-formation. 
We can then find the luminosity for any star-forming galaxy at a fiducial energy of E = 1 
GeV as, 



E L, 



E L 



B,MW 



1> 



^MW 

10 40 GeV 2 sec" 1 GeV" 1 



if, 



1 M yr" 1 



(5) 



where ipMW is the Milky- Way star-formation rate. 

From eqns. ([I]) (j2J), and (jSJ), we estimate the IC contribution to the EGB intensity at 
1 GeV as 



E 2 I E 



k 3x 10~ 8 GeV cm 2 sec" 1 sr" 1 



(6) 



where p+(l) = ~*Mra. i ~ 0.1 M yr 1 Mpc 3 is the cos mic star formation rate at z — 1 
( iHoriuchi et al.ll2009l ). For comparison, the pionic model of lFields et al.l ( 120101 ) gives an inten- 



sity E 2 JelTT-^iGeV ~ 3x 10 7 GeV cm 2 sec 1 sr 1 , with a factor ~ 2 uncertainty in normal- 
i sation. The EGB measured by Fermi -LAT is E 2 Ifijobs l 



GeV 



6 x 10 7 GeV cm 2 sec 1 sr 1 



( lAbdo et al.ll2010j ). Thus, the pionic component alone dominates the overall amplitude from 
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the star-forming galaxies, which is an important contribution to the total observed flux. Not 
evident from our estimate is that the shape of the star-forming spectrum improves upon 
addition of the IC component, as we will see in £JHJ 



3. Targets: Interstellar Photon Fields 




Fig. 1. — The interstellar photon field (ISRF) energy density, shown as a function of the wave- 
length. Th e model include s ther mal components for starlight (optical), dust (IR) and the 
CMB as in lCirelli &: Pancil ( 120091 ) . It includes another thermal component for the starlight 
in the UV. The black (fiducial), red and blue curves represent the curves for (R,h) = (0,5 
kpc), (4 kpc,0) and (0,0) respectively. 



Cosmic-ray electrons in a galaxy will encounter a background photon field, commonly 
called the interstellar radiation field (ISRF). The ISRF is a rich function of both energy 
and geometry. We simplify this complex reality by treating a galaxy as a single zone for 
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the purposes of cosmic-ray propagation and gamma-ray emission. Thus, our ISRF is meant 
to give a sort of effective volume average of the radiation field encountered by cosmic-ray 
electrons. Our choice of ISRF follows those of models for the Milky Way, which we adjust 
at other cosmic epochs via scaling arguments. These scalings use the cosmic star formation 
rate p*{z), whose change with redshift we characterize by the dimensionless function 



S(z) 



P*(0) 



(7) 



where today S(0) = 1. 



For the ISRF at the presen t day (redshift z = 0), we adopt but extend the model pro- 
posed by ICirelli &: Pancil (120091 ) . which consists of three thermal components corresponding 
to the infra-red (or dust), optical, and the Cosmic Microwave Background (CMB). To this 
we add a fourth thermal component for the UV. Our goal is to build a simpl e model that 
can re produce the inverse Compton luminosity of the Milky Way according to IStrong et al. 
(120101 ); their calculation is base d on a far more detailed, spatially-dep endent ISRF as imple- 
mented in the GALPROP code Jstrong et alihoool : Iporter et alJbooa ). The strengths of the 
UV, IR and the optical relative to the CMB decreases with Galactocentric radial distance 
R and also with height h above the Galactic plane. This is seen in Fig. [TJ 



We may express the ISRF specific energy density du lsr i/de = edn lsr t/de as 



A- 



du 



isrf 



du 



isrf 



dX 



de 



7T 2 exp(e/Ti) — 1 



(8) 



a sum of Planck terms, with dimensionless weights f- x for the different background compo- 
nents: % = UV, optical, IR and the CMB. The UV and optical component are from starlight, 
and the IR comes from starlight reprocessed by dust. Our model, as seen in the Fig. [1] 
thus includes the effect of the major components in the GALPROP ISRF, which dominate 
the energy density. For the IR, we tried a non-thermal model that includes effect of both 
scattering and absorption by dust particles, but the effects on our gamma ray luminosity 
model are insignificant. Therefore, we retain the thermal model for IR. We also neglect the 
details at the optical-IR transition, where nonthermal lines appear but whose contribution 
to the energy density is sub-dominant. 

Figure [1] shows our version of the ISRF at different loca tions in the Galaxy. The optical 
and dust components are taken from lCirelli fc Pancil (120091 ). We add another thermal com- 
ponent for the UV. The black, red and blue curves correspond to the positions (R, h) = (0, 5), 
(4,0) and (0,0) kpc. The peaks due to starlight and IR change with position in the Galaxy, 
while the CMB contribution peaked at A ~ 1000 microns remains the same. At the Galac- 
tic centre, the energy densities of these latter components are greater than the CMB, with 
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optical photons dominating. Above the Galactic plane, the CMB becomes comparable in 
energy density and eventually dominates at very large h. 

Table [TJ lists the temperatures of the various ISRF components along with their relative 
weights at zero redshift, /i(0). We chose the relative strengths for regions corresponding 
to 10° < |6| < 20° and 0° < I < 360°, which corresponds to the '10-20' model from 
Cirelli &: Pancil (120091 ) as our fiducial model. This is for the position (R, h) = (0, 5) Kpc 



represented by the black cur ye in Fig. [Tj The U V component at this position is scaled from 
that in the Galactic Centre, (IPorter et al.ll2008l ) in the same proportion as the optical com- 
ponent in these positions. This choice of model, with the overall n ormalisation adjuste d, 
reproduces the Milky Way IC luminosity as a function of energy from IStrong et al.l (120101 ) . 



Table 1: ISRF Parameters 





UV 


Optical 


IR 


CMB 


fiducial /;(0) 
Galactic Centre /i(0) 
T(0)[K] " 


8.4 x KT 17 
1.6 x KT 15 
1.8 x 10 4 


8.9 x KT 13 
1.7 x KT 11 
3.5 x 10 3 


1.3 x KT 5 
7 x KT 5 
41 


1 
1 

2.73 



Because we are interested in inverse Compton emission from star-forming galaxies over 
all of cosmic history, we must specify the cosmic evolution of the ISRF. This will depend on 
which stars c ontribute to starligh t and dust scattered light. The detailed pr esent-day Milky 
Way m odel in lStrong et al.l (120001 ) uses the stellar luminosity functions from lWainscoat et al. 
(Il992h . The UV and IR components derive from short-lived massive stars, and the resulting 
starlight density will simply scale with the star-formation rate. Thus, we adopt the scalings 



h(z) = fw(z) = /uv(0) 3(z) 
f 3 (z) = f m (z) = / IR (0) S(z) 



(9) 



The optical component is less trivial and evolves differently for different stellar popula- 
tions. Massive supernova progenitors, as well as AGB progenitors, are short-lived compared 
to > lGyr timescales for cosmic star formation. Consequently, the optical radiation density 
of these stars would scale as 



h{z) = U(z) = /o P (0) S(z) , (10) 

in step with the UV and IR scaling. For reasonable initial mass functions and star-formation 
histories, these stars should dominate the optical ISRF. But for lower-mass main sequence 
stars and the red giants they become, their cosmologically long lifetimes mean that their 
starlight density at any epoch would scale as the integrated star-formation rate. Therefore, 
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their contribution to the optical component would scale as 



/op(^) = /op(0) J - S , (11) 



This would give a lower optical component and decrease the IC intensity as we will see. We 
will adopt eq. ( fTUl) for our fiducial model, but we will explore the effects of using eq. ( fTTT) ; 
these cases will bracket the true evolution. 

The redshift dependence of the CMB is precisely known and can be expressed entirely 
in terms of its temperature, 

T C mb(z) = (1 + z) T C mb(0) (12) 
with ^(0) = /cmb(0) = 1. The CMB energy density thus grows rapidly as w cm b oc (1 + z) 4 . 



4. Projectiles: Cosmic-Ray Electron Source and Propagation 

An electron with energy E e = 7m scattering off a background photon of energy e will 
produce IC photons with energies around Ej ~ 7 2 e. In order to produce E 1 ~ 100 MeV 
to ~ 300 GeV gamma rays around the Fermi range, electrons with energies of few GeV to 
few TeV are required, depending on the energy of initial un-scattered photon. These are the 
electrons of interest to us. 



High-energy cosmic ray electrons obey a complex transport equation (e.g.. IStrong et al. 



20071 ). This in general would include effects of diffusion, convection, escape, radiative losses, 
etc. However, our approach is to create the simplest one-zone model that captures electron 
source spectrum and their most important losses. 

We model cosmic-ray electron injection via the source spectrum from the Plain Diffusion 



(PD) model of IStrong et al.l (120101 ) . In this model, electrons are accelerated with a broken 
power-law injection spectrum having q e (E e ) oc E" 1 - 8 up to a break at -E e ,break = 4 GeV, 
and above this q e (E e ) oc E~ 2 - 25 . Electrons below -E e ,break can produce inverse Compton 
gamma rays at the low end of the Fermi range, when upscattering optical and UV photons. 
Gamma rays at higher energies are produced by electrons above the break energy. Thus, it 
is important to include the injection break in our calculations^] 



Subsequent to the Strong et al. ( 2010l) analysis of the global Milky Way luminosity, Strong et al. ( 2011 ) 



presented a revised electron injection spectrum that better fits new cosmic-ray electron data. They slightly 
change the power-law indices we have used, introduce a new spectral break at 50 GeV, and an exponential 
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With an interstellar diffusion coefficient, Do ~ 10 28 cm 2 sec _1 (jStrong et al.ll2010l ). the 
diffusion length for electrons over the Compton loss time-scales is i^m ~ 0.01 — 1 kpc, 
depending on the energy. We thus see that electrons will not venture far from the Galaxy, 
but can range vertically outside the Galactic midplane. Hence, very few el ectrons will be 



lost to escape, in strong contrast to ions where escape dominates (|Pohllll994j ). We thus will 
ignore electron escape in our analysis. 

However, since the electrons do propagate throughout and above the stellar and gaseous 
disk, the electron population sees a wide variety of radiation fields. Our single choice is a 
crude approximation meant to represent a sort of average ISRF. As a result, the electron 
propagation equation takes the "thick target" form 

-0 e (£ e ) = q e {E e ) + — [b(E e )<f> e (E e )) (13) 

where q e {E e ) is the source term and b(E e ) = \dE e /dt\ is the total rate of energy loss. The 
equilibrium (i.e., d t (p e = steady state) solution of the cosmic ray flux is 

= fc( * J = W (14) 

As usual in the thick-target limit, the electron flux is directly proportional to the energy- 
integrated injection rate q e (> E e ), and inversely proportional to the total energy loss rate. 
In general, the total energy loss rate of the electrons is given by sum of the inverse Compton, 
synchrotron, bremsstrahlung and ionisation losses, 

b tot (E c ) = b ic (E c ) + b sync (E c ) + b hlcm (E e ) ; (15) 

we now address each in turn. The ionisation loss is unimportant throughout the entire 
electron energy range of interest here. 

The inverse Compton loss in the Thompson regime is given by the simple expression 

, Thomps 4 ( E e x - 

b ic (E e ) = -a T c U- mT f - 

3 \m c 2 1 

« 1 x 10- 8 GeV sec- 1 ( ^-^] ( — ) (16) 

V 1 eV cm" 3 J V 10 TeV ' 



cutoff at 2 TeV. We have verified that, when keeping fixed the injection at the 4 GeV break, the lStrong et al. 
(|201ll ) spectrum does not strongly affect the overall IC amplitude, and over the Fermi EGB energy range 
the resulting IC luminosity is similar near 200 MeV, and slightly lower at higher energies. 
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However, the full Compton cross section takes the Klein-Nishina form described in more 
detail in §5.11 Corrections to the Thompson limit become important when A e = 4 e'j/m c 2 3> 
1; in a starlight-dominated ISRF with e ~ 1 eV, this occurs for electron energies 7 3> 10 5 . 
This energy regime is important for our calculation and therefore, we calculate the inverse 
Compton loss rate by combining the Klein-Nishina cross section self-consistently with the 
(redshift-dependent) ISRF: 

, far 

\ J - J erti 



b ic (E e 



de — - — e 



de 



dE P 



dE P 



3 0"t c 



dn[ sr f 

de ^T (e) 



(4 (^) 2 -Ae)g-l 



(m c 2 /2E c )' 2 



x 



2q ln(g) + (1 + 2q)(l - q) + " q) 



(17) 



as m 



Cirelli fc Pancil (120091 ) . where 



E P 



4e7 2 (l - E em /~/m c 2 



(Felten 


1965; 


Pohl 


1993) 



Compton losses are all important, depending on the part of the spectrum. For electron 
energies lower than ~ 1 GeV, bremsstrahlung is dominant, with Compton losses being 
comparable close to the break energy at E c = 4 GeV. Therefore, instead of seeing two 
breaks corresponding to these two cosmic ray propagation regimes, only one is seen near the 
break energy as described in §H 

Note also that for the highest-energy electrons, inverse Compton losses b ecome catas- 
trophic, i.e., in e ach scattering event the fra ctional changes in the electron energy approach 



5E e /E e ~ 1 (e.g.. lBlumenthal fc Gouldlll970l ). In this case the losses can no longer be treated 
as a smooth function bi C (E e ) that averages over many scatterings each with 5E e /E e 1; we 
do not include these effects, which impact gamma-ray energies of E cm > 1 TeV. 



The other energy loss rates are given as follows (IHayakawal Il973t lGinzburglll979l ): syn- 
chrotron losses 

2 1 / r-i \ 2 
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m cr 
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o 7r \m 



3 x 10" 10 GeV sec" 1 
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E P 



(19) 



are controlled by the interstellar magnetic energy density U mag . The bremsstrahlung losses 
depend on whether the medium is ionised or neutral. Due to the presence of both H I and H 
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II, the bremsstr ahlung losses are due to both neutral and ionised hydrogen (lHayakawalll973 
GinzburejE979[ ): 
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(23) 
(24) 



Here nn i and % n are the number density of relativistic electrons in the interstellar medium, 
which is equal to the number density of H I and H II. 

Different loss mechanisms dominate at different energies. Inverse Compton losses scale 
as b lc oc E 2 e in the Thompson regime, but for large E e this drops off to a logarithmic energy 
dependence due to Klein-Nishina suppression. Synchrotron losses have 6 sync oc E 2 , and thus 
are proportional to and are comparable to inverse Compton losses at moderate energies 
E e ~ lGeV — TeV. At high electron energies, inverse Compton is suppressed and the 
synchrotron losses dominate. Below E e ~ 1 GeV, bremsstrahlung losses &brcm oc E e become 
important. 

The losses themselves evolve with redshift, through the dependences on interstellar den- 
sities, i.e., the ISRF energy density U isT {, the interstellar magnetic energy density ?7 mag = 
B 2 /8tt, and the number densities nj sm of interstellar particles. The magnetic field couples to 
all cosmic rays, for which the ions dominate the energy density. This coupling is likely re- 
sponsible for the approximate equipartition of magnetic field and cosmic-ray energy densities 
in the local interstellar medium. Therefore, we assume that the magnetic field energy density 
scales with the cosmic-ray ion flux, and therefore with the star-fo rmation rate: Um^ oc S(z). 
For our fiducial model, we use B = 4 /i G, to match the mo del of IStrong et al.l (120101 ). which 
is comparable with the typical value in IStrong et al.l (120111 ) . 



We take interstellar particle densities, which control bremsstrahlung losses, to scale as 



(25) 
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This c an be viewed a consequence of the disk radius scaling -Rdisk oc (1+z) in the lFields et al. 



(120 lOl ) model, and to reflect proportionality between galaxy disk and dark halo sizes. The 
fiducial number densities of neutral and ionised hydrogen, i.e. n H 11,0 — i,o = 0.06 cm -3 . 
The helium content of ionised and neutral gas is also included, with 7/ He = 0.08. These 
parameters are broadly consistent with values at intermediate heights h above the Galactic 
plane, as befits the volume average represented by our one-zone model. For reference, the 
values of n ism and B used at the Galactic centre are much higher at B = 8.3 /iG and 
n H n,o = ™H i,o = 0.12 cm" 3 . 



5. Inverse Compton Emission from Individual Star-Forming Galaxies 

As seen in eq. (j2J), we first compute the luminosity of single galaxy and then average 
suitably over the luminosity function to get the luminosity density or cosmic volume emis- 
sivity. In this section, we address the former, focussing first on the IC emissivity within the 
volume of a single galaxy, then using this to compute the total luminosity from that galaxy. 



5.1. Inverse Compton Emissivity of a Star-Forming Galaxy 

The specific IC volume emissivity within a galaxy is the rate of producing gamma ray 
photons, o?r e - 7 _ s . e - 7 /(ii? em per background photon, multiplied by the ISRF number density, 

"isrf 

d(jry dN^ y cTg— 7 _^ e - 7 dn 1ST f (26) 



dE em dVdE em dt J dE era de 

Here, the gamma ray rate per unit ISRF photon is a product of the cosmic ray flux and the 
IC cross section 

d^e~ ■y— ►e - 7 d<7 e -^,— > e — -y r n i-r\ 

~^TF = ~^VF ^ ■ ( 2? ) 

uc em u, - C/ em 

that in general takes the Klein-Nishina form. The full Klein-Nishina cross section is a com- 
plicated function of the energies involved. Moreover, the inverse Compton process produces 
anisotropic emission if either the cosmic-ray or ISRF populations depart from isotropy. While 
the cosmic-ray electrons of interest are well-approximated as isotropic, in a real Galaxy the 
photon field is certainly anisotropic as wel l as spatially-varying. These effects are included 



in GALPROP (IMoskalenko fc Strong! |2000| ). but we will ignore them in our simple one-zone 



galaxy model. This is a ~ 20% effect which will not dominate our final error budget. 
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For an electron with Lorentz factor 7 = E e /m e , the Thompson limit is given as 

4 e7 < mc 2 (28) 

For the electron spectra we consider, the Thompson limit is a good approximation for CMB 
and IR photons, but fails for UV and som e op tical photons. Hence, we use the Klein-Nishina 



result via the prescription of Ijonesl ( 119681 ) and iBlumenthal fc Gouldl (Il970l ). This differential 
cross section is expressed as, 

^ (6, E.,q) = l2L( 2q ln(q) + (1 + 2q)(l - q) + - q) 



dE em y ' 4e 7 2 V 21 + A e g 

where a? is the total Thompson scattering cross section and q appears in eq. (fT8"j) . In the 
Thompson limit, the last term in the above equation become negligible. 

The inverse Compton volume emissivity within a one- zone galaxy is thus 

d qi _ f & dr - ^ e - 7 dn isrf = r d dn*a r dEe< p e ( Ee) da *-^-~< (29) 



dE cm J dE em de J de J dE L 

We evaluate the emissivity numerically, by using the ISRF, cosmic ray flux density, and the 
cross section mentioned above, §3] and §HJ That ISRF density n isr f is the same as in the 
inverse Compton energy loss equation eq. (EEL so that we calculate the two self-consistently. 



It will be useful to cast the emissivity in the form 



f A r tfajsrf , f 1 dT(t,E am _) \ 

dq-f TT J ae de e {e dE cm ) _ / 1 dT(e, E em ) \ / 1 q c (> E e ) da x 



dE em lsrt Jde^jfe " \e dE cm / mTt \e b dE, 

where the angle brackets indicate an average weighted by the ISRF energy density distribu- 
tion, and where the differential and total ISRF energy densities are w isr f = e dn isr{ /de and 
^isrf = / de Ui S rf respectively. We see that the amplitude of the emissivity scales as 

^isrf ^ic 
9 7 ^~1e ~ y^e 

Therefore, the photon emissivity, i.e., the gamma-ray source rate, is proportional to the 
electron source rate, but also depends on an appropriately weighted ratio of the inverse- 
Compton losses to the total losses. This ratio controls the shape of the spectrum. If the 
losses are Compton-dominated, then b ~ b lc and thus, the two source rates amplitudes are 
directly proportional, g 7 oc q e , and the IC emissivity is independent of the ISRF energy 
density, U- mT f. Thus, the number and energy of inverse Compton photons is proportional to 
that of cosmic-ray electrons, because in a steady state, the energy injected into electrons 
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must equal the energy emitted in IC photons. This situation defines perfect calorimetry. On 
the other hand, the ISRF shape in energy space is what weights the average in eq. (130]) . and 
thus the shape of the gamma-ray spectrum does remain sensitive to the ISRF even in the 
case of perfect calorimetry. 

However, in reality, synchrotron and bremsstrahlung compete with inverse Compton 
losses at different energies, as seen in §H In the Thompson regime, the ratio of synchrotron- 
to-inverse Compton losses 

frsync Thomps Ub_ _ „ . ( 1-1 eV Cm' 3 \ / B \ 2 
&ic C/isrf ~ V ^isrf J WSJ ' 1 J 

is energy-independent. For our fiducial magnetic field value, this ratio is not far from unity 
and hence, the shape does not change much for different parameters describing the ISRF 
and the magnetic field. On the other hand, for E e > 1 TeV, the IC losses are in the Klein- 
Nishina limit and grow only logarithmically in energy; then the synchrotron losses dominate 
In this case, we have q 1 ~ (U- mT {/U m3 ,g)q e . Thus, for high-energy cosmic-ray electrons, the 
calorimetric approximation increasingly breaks down. 

The competition among losses depends on the way the interstellar densities evolve with 
redshift. Short-lived, massive stars dominate most of the ISRF except possibly for some of 
the optical range, and so the energy density scales with the star- formation rate Ui ST { cx S(z). 
Due to cosmic-ray ion equipartition, we take t/ mag oc S(z) as well, and thus the Ui sr f/U mag 
ratio doesn't change dramatically throughout the cosmic history. Similarly, we also include a 
strong evolution n ism oc (1 + z) 3 evolution to the interstellar particle density, which increases 
rapidly towards z ~ 1, roughly in step with the cosmic star-formation rate. Thus, the 
inverse-Compton/bremsstrahlung ratio of energy losses remains roughly constant as well out 
to z ~ 1. The net result is that the spectral shape of inverse Compton emission does not 
evolve dramatically with redshift in our model. 



5.2. Total Inverse-Compton Luminosity from a Single Galaxy 



As seen in knowing the specific inverse Compton emissivity dq^j dE em within a single 
star-forming galaxy, we can calculate the specific IC luminosity of the galaxy as 



L 7 (E e 



dNjtc 
dt dE om 



dq 1 
dEL 



-dVi 



ISM 



(32) 



and we have seen that g 7 oc (bi C /b to t)qe- Following Fields et al.l (120101 ) . we take supernovae 
as the engines of cosmic-ray acceleration. Thus we scale the electron injection rate with 
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Fig. 2. — The inverse-Compton gamma-ray luminosity spectrum for a single galaxy with 
the Milky- Way star-formation rate. The solid black curve represents our fiducial model, 
with an ISM density n H 1,0 = ii,o = 0.06 cm -3 , magnetic field, B = 4 \iG. This is for 
the ISRF consistent with (R, h) = (0, 5) kpc. This model is des igned to provide a good 
fit to the GALPROP Plain Diffusion model of IStrong et al.1 (120 lOl . dot-dashed blue curve). 
The effect of some variations of gas density, magnetic field and ISRF are shown in red. 
The solid red curve shows the effect of reducing B to 2 \i G. The dotted red curve is for 
n H i,o — n H n,o — 0.03 cm -3 . The red dashed is a toy model with the ISRF including the 
CMB alone. The solid blue shows the single galaxy template for the Galactic Centre, with 
nn i,o = "^h n,o = 0.12 cm -3 , and B = 8.3 fiG. 
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the supernova rate, i.e., we have L 7 ~ J q e dVisM oc -Rsn oc ip. This linear proportionality 
captures the main dependence of IC emission on star-formation rate. 

Of course, in reality the cosmic ray flux and the ISRF are both functions of position 
in the galaxy. But we are describing a one-zone model, in which our choices of the fiducial 
p arameters are made in order to provide the best match to the GALPROP Milky Way model 



m lStrong et all (120101) . Specifically, w e have created a best-fit model that mimics well the IC 
luminosity in the IStrong et all (120101 ) Plain Diffusion case. We adopt the same cosmic ray 
electron injection spectrum and amplitude, which are themselves chosen to provide a good 
fit to the observed local cosmic-ray electron spectrum. The interstellar parameters are the 
ISRF component amplitude fi, the interstellar particle density n ism , and the magnetic field 
strength B\ physically, these representing an average over the volume occupied by cosmic-ray 
electrons. 

Thus, we have a single-galaxy IC luminosity that takes the separable form 

L 7 (ip,z) = -^—L MW (z) (33) 

WW 

where we have suppressed the energy dependence for clarity. Here we explicity display the 
dominant sensitivity to star-formation, via the overall linear proportionality to the star- 
formation rate. The Milky- Way model Lmw( z ) is f° r a galaxy with the present-day Milky- 
Way star-formation rate, but allowing for redshift evolution of the interstellar matter and 
radiation densities. 



Our IC luminosity models for z = are shown in Fig. [2J along with that of IStrong et al. 



(120101 ) Our fiducial model described in $3j §4] and §5. H is represented by the solid black curve. 
The main feature of the inverse Compton E^ m Luw spectrum is that it has a broad maximum 
corresponding to the break in the electron injection spectrum at -E e ,break = 4 GeV. In our the 
fiducial model, the optical component dominates the ISRF, and thus we expect the injection 
peak to correspond to an IC peak at E em ~ 7 2 e ~ 100 MeV. This is about what our fiducial 
model predicts and is consistent with GALPROP. 



Our fiducial model does a good job of reproducing the IC results of IStrong et al.l ( 120101 ). 
over most of the Fermi energy range. Thus, our results are equivalent to normalizing to the 
GALPROP luminosity for galaxies with the Milky Way's star formation rate. We do not 
reproduce GALPROP's steep decrease in the spectrum beyond 100 GeV, which may result 
from our neglect of the catastrophic nature of the electron energy losses in this regime. A 
fitting function for our fiducial IC model appears below in Table [2J the fit is good to better 
than 2% over the Fermi energy range. From this fit we see that at 1 GeV the logarithmic 
slope is around —0.1, so that Lmw ~ Ecm' 1 - This is considerably flatter than both the pionic 
spectrum and the observed EGB at these energies. 
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Our choices of interstellar parameters for the fiducial model were made to provide a good 
fit to the GALPROP luminosity, and also are physically reasonable as volume averages. 
But of course other parameter choices are possible. To illustrate the sensitivity to these 
parameters, Fig. [2] shows variations of the luminosity with interstellar particle density nj sm 
and magnetic field strength B. Decreasing the gas density suppresses the bremsstrahlung 
losses relative to the total. This increases the ratio of b\ c /b\, rerai which governs the IC spectral 
shape at lower energies. The result is a boost in the IC output at lower energies, which also 
shifts the IC peak towards lower energies. At high gamma-ray energies beyond ~ few GeV, 
synchrotron losses dominate, and thus the ratio bi C /b syQC controls the spectrum; reducing the 
magnetic field increases the luminosity. 

The choice of ISRF is also important. Our fiducial model adopts an ISRF consistent 
with the GALPRO P (R, h) = (0, 5) kpc emission, as it best reproduces the IC luminosity of 



(IStrong et al.ll2010l ). If we change the ISRF to that of the Galactic centre, the IC signal is 
amplified, and the peak shifts to lower energies, as shown by the solid blue curve in Fig. [2j 
At the Galactic centre, the magnetic energy density is increased by a factor of 4 relative to 
the fiducial case, but the photon energy density goes up by more than a factor of 10. So both 
the fractions, b ic /b and b sync /b are higher, but so is b ic /b sync . And so the peak occurs at a 
lower energy, but has a greater amplitude. Finally, for illustration we turn off all components 
of the ISRF except the CMB, which substantially reduces the ISRF energy density and thus 
degrades calorimetry. We thus find that the overall signal is reduced and the peak shifts to 
far lower energies, as expected; indeed this peak is now off scale in Fig. [2J Of course, such a 
model is unphysical, because we are studying electrons born in star-forming systems where 
starlight is always present by definition. 

Thus, there are several uncertainties in our simple model of the gamma ray spectrum 
for a single, Milky- Way like galaxy, such as choice of the average ISRF, the electron density, 
and the interstellar magnetic field. However, for reasonable interstellar models, these all 
lead to modest changes in the normalisation of the luminosity, at most tens of percent. 
Larger chang e s wou ld be possible if one were to depart from the single galaxy template in 



Strong et al.l (120101 ). which we try to reproduce and that is based on a wealth of Milky- 
Way data. Allowing oneself sufficient freedom, the inverse Compton EGB signal can be 
adjusted by changes, e.g., in the U^/IIb and U isr f/n ism ratios that control the degree of 
electron calorimetry. Of course, this would then drive the system away from the rough 
energy equipartition observed in the Milky Way. Even then, large increases in these ratios 
would only increase the completeness of electron calorimetry and would raise the IC signal 
by a factor < 2; on the other hand, large decreases in these ratios would spoil calorimetry 
and could substantially lower the IC luminosity. Thus, the cosmological prediction of the IC 
contribution of the star- forming galaxies to the EGB should be fairly robust, unless there 
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are departures from the equipartition implicit in the Milky- Way-based normalization. 



6. Results: Inverse Compton Contribution to the Extragalactic Background 



Having established a one- zone galaxy template, and explored physically plausible varia- 
tions in it, we proceed to compute the IC contribution to the EGB measured by Fermi -LAT 
from star-forming galaxies. This is given by the well-known expression 
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47T 
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(34) 
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Here we assume a flat ACDM universe 



with H = 71 km s _1 Mpc -1 , ft 



matter 



0.3, and Q A = 0.7, following Fields et all ( 120101 ). 



The cosmic IC luminosity density (i.e,. cosmic volume emissivity) is given by substitut- 
ing our single-galaxy luminosity of eq. ( 133|) into the luminosity function of eq. ([2]). Because 
we have L lc oc if), the IC emission of a star-forming galaxy traces its star-formation rate. 
Thus, the IC luminosity function dn za \/dL^ at each redshift is proportional to the lumi- 
nosity function of any tracer of star formation at that redshift. Suppressing the energy 
dependence for clarity, we have 
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(35) 
(36) 



where the P * = (j^ ^gai) = / j> dn „^ / dL^ dL^ is t he cos mic star-formation rate. Following 



Fields et al.l (120101 ) we adopt the iHoriuchi et al.l (120091) cosmic star-formation rate, and a 
Milky Way star-formation rate ^mw = l.O7M /yr (IRobitaille fc Whitneyll2010l ). 



The overall amplitude of the cosmic IC luminosity density is thus linearly proportional 
to the cosmic star-formation rate. Of course, the detailed spectral shape of the luminos- 
ity density reflects redshift dependence of the interstellar matter and radiation encoded in 
Lmw( z )- 

Figure [3] shows the IC contribution to the EGB for several choices of the Milky- Way 
spectrum Lmw- The results for the fiducial model appear as the solid black curve, and 
incorporate the full redshift-dependence of interstellar density encoded in L MW (z). A fitting 
function for the fiducial model appears in Table [2J valid over the Fermi energy range. For 
comparison, Fig. [3] shows the pionic contribution in red. The IC contribution is about 10% 
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Fig. 3. — The inverse-Compton EGB intensity from star- forming galaxies as a function of 
energy. The solid black curve represents our fiducial IC model, using a single-galaxy model 
where interstellar densities evolve with redshift. The dotted and the dashed curves show the 
single galaxy template that does not evolve but is fixed to the behavior at z — 0, 1 and 2. 
The solid blue curve is for an evolving single galaxy template with an ISRF corresponding 
to the Galactic centre model. The dotted blue shows the effect of scaling the optical ISRF 
with the cosmic stellar mass dens ity feci. [TTj). The r ed curve shows the pionic signal from 
normal star- forming galaxies as in iFields et al.l ( 2010 ). 
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of the pionic contribution at the ~ 300 MeV peak in E 2 1\ growing to about 20% at 10 GeV. 
This is consistent with our expectations from the order-of magnitude calculations of £j2j 

The shape of the IC curve is rather smooth, and rather flat in E 2 I. This is because the 
spectrum it is a redshift-smeared version of the single-galaxy IC emission. This stands in 
contrast to the pionic curve that displays its characteristic peak and a steep ~ E~ Sp dropoff 
at large energies, with s p = 2.75 the propagated proton spectral index. 

We have verified that the bulk of the IC signal comes from reshifts around z ~ 1. That 
is, the integrand in eq. ( jMJ) peaks in this regime. A comparison of the corresponding curves 
in Figs. |2] and [3] shows that the IC peak in the galaxy luminosity rest-frame energy E em ~ 100 
MeV is transformed into an EGB IC peak at ~ 40 MeV. This translation in energy space is 
consistent with redshifting by a factor 1 + z . 

The sensitivity to the interstellar density evolution adopted in our fiducial model is 
explored by the dashed and dotted black curves in Fig. [31 These show the effect of assuming 
a single galaxy spectral template that does not evolve with redshift, but rather at all epochs 
uses non-evolving interstellar densities corresponding to our model Lmw(^o) evaluated at 
some fixed redshift zq. Results are shown for this non-evolving source spectrum fixed for 
the entire cosmic history to the spectrum for redshifts Zq — 0, 1 and 2. It is clear that the 
actual cosmological result shown by the solid black curve is closest to the zq = 1 case, where 
the bulk of the IC signal originates. Moreover, we see that all of the single-redshift source 
spectra lead to EGB spectra very similar to the fiducial evolving Lmw{ z )- This reflects the 
insensitivity of our source spectrum with redshift, which itself arises due to calorimetry and 
the near-constancy of the f/; sr f : U mag : n density ratios. 

To give a sense of the effect of the adopted ISRF, we have computed the IC EGB 
resulting from an evolving ISRF whose spectral shape is appropriate for the Galactic centre. 
This model corresponds to the solid blue curve in Fig. [2j Here again, the main result is 
that the curves are very similar. In detail, the Galactic centre has a slightly higher ratio 
Uisrf/U mag of starlight-to- magnetic energy density. Thus, electron energy losses are more IC 
dominated, and finally the overall IC flux is closer to calorimetric and so slightly higher. 

Our fiducial model assumes the ISRF components (other than the CMB) are dominated 
by short-lived stars and thus have a redshift evolution that scales with the cosmic star- 
formation rate. An extreme alternative is that the optical ISRF is dominated by long-lived 
main sequence stars whose energy density scales with the stellar mass and thus the integrated 
cosmic star- formation rate, as in eq. (ITTj) . The dotted curve in Fig. [3] shows the EGB for 
this case. The signal is somewhat smaller than the fiducial model across all energies. This 
is easily understood: at high redshifts the star-formation rate is higher than today, but the 
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Fig. 4. — A summary of the star-forming galaxy contributions to the EGB inte nsity, shown 



as a f unction of energy. The blue points represent the Fermi -LAT data points ( lAbdo et al. 



2010|) . The solid red and blue curves show the gamma ray intensity due to pionic and fiducial 



IC components from star-forming galaxies as in Fig. [3j The solid black curve represents the 
sum of the two. The shaded band gives our estimate of the uncertainty in the predicted 
total signal, which is dominated by systematic errors that are common to the IC and pionic 
components. 
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stellar mass is lower and thus offers less optical photons and reduced calorimetry. But the 
effect is not dramatic; even in this extreme case, the reduction in the IC signal is less than 
a factor of two relative to our fiducial model. 

Thus, we see that the IC results do depend on the adopted interstellar densities and their 
redshift evolution. But for reasonable choices, the variations in the final IC EGB result are 
small. The IC calculation in this sense appears rather robust to the systematic uncertainties 
in the modeling. 



It is useful to compare the IC luminosity density in eq. fl36|) with that for pionic emission. 
As noted above, the cosmic IC emissivity is proportional to the cosmic star-formation rate, 
alo ng with a weak sen sitivity to the redshift evolution of interstellar densities. In contrast, 



the Fields et al.l ( 120101 ) pionic emission depends on the product of cosmic star-formation rate 
and mean interstellar gas mass, via the nonlinear scaling L n ~ ip 1+ul with u = 0.714. Thus, 
the pionic emissivity represents a different moment of the cosmic star-formation rate C n = 
(L^rigai) ~ (ip 1+UJ n ga \). This nonlinear moment gives different results depending on how the 
star-formation distribution evolves with redshift. Two extreme limits correspond to: (a) pure 
density evolution, wherein galaxies have constant star-formation rates but evolving number 
density; versus (b) pure luminosity evolution in which the star-forming galaxy density is 
constant but the star-for mation rates all evo lve with redshift. The pionic curves we show 



correspond to the fiducial iFields et al.l ( 120101 ) pure luminosity evolution case, which gives a 



contribution to the EGB larger by a factor ~ 3 than in the pure density evolution case. In 
the IC calculation, the pure luminosity and pure density evolution are degenerate, in that 
these models all give the same cosmic star- formation rate and thus the same IC output. 

Figure H] represents the central result of this paper, summarizing the EGB predictions 
for gamma-ray emission from star-forming galaxies in the Fermi energy rang e. We show 



the fid ucial IC curve from Fig. ([3]), the fiducial normal-galaxy pionic model of IFields et al. 



(120 101 ). and the total EGB intensity that sums these components. Over the 100 MeV to 
300 GeV range, the pionic component dominates the normal galaxy contribution. But as 
we have seen, the IC curve is a much flatter function of energy. The IC component thus 
becomes systematically more important at higher energies. Consequently, high-energy slope 
of the total emission is less steep than that of the pionic signal, and closer to the Fermi-LAT 



data shown in blue (lAbdo et al.ll2010fl ). 



Including the IC emission thus improves the agreement between sta r-forming prediction 



and the observed EGB slope, which had been a weaker point of the IFields et al.l (120101 ) 
pionic-only model. Up to about ~ 10 GeV, the central values of the model fall below the 
data, but are consistent within the theoretical and observational error budget (discussed 
below). At higher energies, the model underpredicts the data. Our pionic model neglects 
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starb urst galaxies; it m ay well be that their pionic emission is important at this energy range 



Lacki et al.l 120111) . In a ddition, unresolved blazars are also guaranteed to play a role 



Stecker fc Ventersll201ll ). In any case, star-forming galaxies clearly are an important 



(e.g., 
(e.g., 

contribution to the Fermi EGB signal, and could well be the dominant component. 

Note also that the IC curve does no t beco me important at low energies until far below 
the pion bump. Thus, the Fields et al.l (120101 ) conclusions stand: the pion feature should 
still remain as a distinguishing characteristic of a significant star-forming contribution to the 
EGB. Measurements below ~ 200 MeV should show a break in the EGB slope if star-forming 
galaxies play an important role. Such a spectral feature provides one way to discriminate 
between star-forming galaxies and other sources, such as the guaranteed contribution from 
unresolved blazars. 



Table 2: Fitting Functions 
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Our estimate of the uncertainty in the total normal galaxy emission is represented by 
the s haded band i n Figu re |H The errors in the pionic contribution dominate, and are taken 
from iFields et al.l ( 120101 ). For the IC contribution, the errors are dominated by systematic 
uncertainties in the Milky- Way star formation rate ^mw and in the normalization p*(0) of 
the cosmic star-formation rate. These factors are common to the overall normalizations of 
both the pionic and IC signals, cf. eq. (|36|) . Hence, the addition of the IC signal to the total 
does no t substantial l y incr ease the error budget, and so we have retained the uncertainty 
band of IFields et al.l (120101 ) . The errors as shown are slightly underestimated at the highest 
energies, but we will see that other effects enter at this regime. 

We have ignored the effect of intergalactic absorption of the high-energy gamma rays via 
77 P hi — > e + e~ photo-pair production in coll i sions with extragal actic background light (e.g., 



Salamon fc Steckerlll998l : lAbdo et al.ll2010d : IStecker et al.ll2012l ). This attenuation starts to 
become significant for gamma rays over few tens of GeV. Therefore, beyond these energies 
our results will be suppressed by amounts that depend on the optical depth for these high- 
energy gamma rays. In this context it is worth noting that Fermi -LAT observations of the 
z = 0.9 active galaxy, 4C +55.17 do not show significant absorption of gamma rays up to 
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about 100 GeV (jMcConville et al.ll201ll ). This object lies within the z ~ 1 regime where 
most of the IC contribution to the EGB occurs; thus we might expect that attenuation 
should be mild at Fermi energies. 



7. Discussion and Conclusions 

We have calculated the contribution to the extragalactic gamma-ray background due 
to inverse-Compton emission from star-forming galaxies. To do this we model the IC emis- 
sion in individual star-forming galaxies, which arises from cosmic-ray electron interactions 
with the ISR F. We normalize to the present-day Milky Way IC luminosity as determined by 



GALPROP ( iStrong et al.ll2010l ). Indeed, we hope our simplified models helps to illuminate 
some of the rich physics in the GALPROP model. We also provide a prescription for red- 
shift evolution of interstellar matter and energy densities, based on equipartition arguments. 
Assuming cosmic rays are accelerated by supernovae implies that a galaxy's IC luminosity 
scales ip; this further implies that the cosmic IC luminosity density or volume 

emissivity is proportional to the cosmic star-formation rate: £ 7 oc p*. 

This linear dependence of IC luminosity with star-formation rate has important im- 
plications in light of the star-forming galaxies resolved by Fermi. Namely, the L[ C oc ip 
trend provides a po or fit for Fermi gal axies, for which the observed correlation is non-linear: 



L ic oc ^ L4±0 - 3 (e.g.. lAbdo et al.ll2010dl ). This implies that the IC component is subdominant 
in Fermi galaxies, adding indirect evidence for the primacy of pionic emission in star-forming 
galaxies and prefiguring the trends we find for the EGB. 

Turning to the EGB, we find that the IC contribution has a very broad maximum in 
E 2 I[ C at E ~ 40 MeV, falling off very gradually away from this peak. This shape is a 
redshift-smeared reflection of the underlying Milky- Way-like spectrum from the individual 
galaxies. In fact, the IC emission is so broadly peaked that is it effectively featureless. This 
contrasts markedly with the distinctive pionic feature from cosmic-ray hadronic interactions 
in star-forming galaxies. 

The amplitude and shape of the IC contribution to the EGB depends on the nature 
of the interstellar radiation and matter fields, and their evolution with redshift. However, 
we find that when we adopt different but physically motivated variations to our fiducial 
model, the final EGB predictions change very little. This rough model-independence of the 
IC EGB is a consequence of partial calorimetry, i.e., the fact that IC photons represent 
a substantial fraction of cosmic-ray electron energy loss, so that the IC output is tied to 
cosmic-ray production via energy conservation. 



-27- 



We find that in all of our models, the IC signal is always smaller than the pionic 
contribution from normal galaxies. This largely follows from our normalization to the Milky 
Way emission, where the GALPROP model finds (and Fermi observations imply) that pionic 
emission dominates over IC. However, while the EGB IC component is smaller, it also has a 
substantially flatter spectrum, so that the IC becomes increasingly important away from the 
pionic peak. This implies that, at least for normal star- forming galaxies, the IC contribution 
should dominate over pionic at high and low energies. However, at higher energies the 
opacity of the universe becomes important and losses become catastrophic, and requires 
different techniques to handle correctly. At energies below the Fermi-LAT range, the IC 
emission will be supplemented by processes such as brems strahlung, which w e have not 



included; a detailed treatment of the MeV regime appears in lLacki et al.l (120121 ) . 



The approximate calorimetric relationship between IC photons and cosmic-ray electrons 
has important consequences. One is that as long as cosmic-ray electrons are accelerated with 
similar spectra and efficiency everywhere, their resulting IC emission will not depend strongly 
on their environment. For this reason, we expect that the IC output per supernova should 
be the same for normal galaxies and starburst galaxies, at least to zeroth order. This would 
mean that, unlike the pionic case, the IC emission from all star-forming galaxies can be 
treated on the same footing. Thus, we have not excluded starbursts in our analysis, as they 
were from the pionic signal. If we did so it would only reduce the overall IC signal from 
normal galaxies to even less than the pionic contribution. 

Another consequence of partial calorimetry is that the main redshift dependence of 
the cosmic IC emissivity is a linear dependence on the cosmic star-formation rate. Thus, 
our results are independent of whether the cosmic star-formation rate is a result of pure 
luminosity evolution, pure density evolution, or something in between. This is in contrast 
to the pionic case, which depends nonlinearly on the star-formation luminosity function and 
so breaks the degeneracy between the pure luminosity and pure density evolution cases. 

To simplify the discussion, the IC and pionic models presented here neglected the effects 
of Type la supernovae, impli citly assuming i nstead that all supernovae in star-forming galax- 



ies are due to core-collapse. iLien &: Fields! (120121 ) considered in detail the effect of Type la 



supernovae on the pionic signal. They found that a self-consistent treatment includes both 
the addition of Type la explosions as cosmic-ray accelerators, but also as part of the Milky- 
Way normalization of the cosmic-ray /supernova ratio. The effects largely cancel, so that in 
the end, the net pionic galactic luminosity and EGB does not change substantially. In the 
case of IC emission a similar cancellation will occur. The only new contribution of possible 
importance comes from Type la explosions arising from long-lived progenitors in quiescent 
(i.e., not actively star-forming) galaxies such as ellipticals; the pionic emission from these 
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systems could be large if they have a substantial reservoir of hot, X-ray-emitting gas. But 
supernova rate in these galaxies represents a subdominant fraction of cosmic Type la ac- 
tivity, which itself is substantially smaller than the cosmic core-collapse rate. Thus, the IC 
contribution from these systems will be small. And so inclusion of Type la supernovae in 
a self-consistent way would change our results very little, less than the uncertainties in the 
model. 

To our knowledge, this paper present the first discussion of the IC contribution from star- 
forming galaxies to the EGBg There remains room to improve on our model. Future work 
would benefit from observational progress in clearly identifying an IC signal from individual 
star-forming galaxies, at energies away from the pionic peak. Theoretical work would benefit 
from a more detailed model for the ISRF and its evolution, and from the use of additional 
multiwavelength constraints on the cosmic ray electrons. 

More broadly, a solid identification and quantification of the main components of the 
EGB remains a top priority for gamma-ray astrophysics and particle cosmology. Extending 
the Fermi EGB energy spectrum to both higher and lower energies will provide important 
new constraints. At sufficiently high energies, the cosmic opacity due to photo-pair produc- 
tion must become apparent if the EGB is dominated by any sources at cosmological distances 
(e.g., ISalamon & Steckerlll9 98; Ab do et al.ll2010d ; IStecker et al.ll2012[ ). And at energies just 



below those reported in lAbdo et al. 



(l2010if ). a break in the EGB spectrum is an unavoidable 



prediction if the signal is dominantly unresolved pionic emission from star-forming galaxies 



(bo th normal and starbursts 



ies (lAndo fc Komatsu 



2012 



2006 



An independent prob e of EGB origin lies in anisotropy stud- 



Ando fc Pavlidoul 120091 ; iHensley et al.l l2010t lAckermann et al. 



Malyshev fc Hoggll201lh . Regardless of the outcome, an assay of the EGB components 



will provide important new information (and perhaps some surprises!) about the high-energy 



cosmos. 
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2 As we were in the final stage of preparing this paper, we became aware of the work of 
which addresses similar issues. 



Lacki et al. 



torn 
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